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Abstract 

The momentum dependence of the empirical scalar and vector potentials 
needed for describing relativistic heavy-ion collisions at 1 GeV/nucleon is 
compared with that derived from self-consistent Dirac-Brueckner calculations 
using the Bonn-A potential. Our calculated scalar and vector potentials ex- 
hibit a weak momentum dependence and their magnitudes decrease with nu- 
cleon momentum, in close similarity with the momentum dependence required 
by the empirical potentials. The effects of explicit momentum dependence on 
the properties of equilibrium nuclear matter are found to be small. 
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Heavy-ion collisions at various energies open new frontiers and provide opportunities for 
nuclear physics. Theoretical description of heavy-ion collisions is usually based on transport 
models. For heavy-ion collisions at beam energies from a few hundred MeV/nucleon to a few 
GeV/nucleon, the reaction dynamics is mainly governed by the stochastic two-body scatter- 
ing and propagation in the mean-field potential. Transport models that include both effects, 
such as Boltzmann-Uehling-Uhlenbeck (BUU) [||], quantum molecular dynamics (QMD) 0, 
and the relativistic extensions, RBUU and RQMD have been extensively used for 
the description of heavy-ion collisions in this energy region. 

In the BUU or QMD calculations, the nucleon potential is usually parameterized in terms 
of Skyrme-type effective interaction The importance of the momentum dependence 

of the nucleon potential was realized when analysing nucleon flow and pion production in 
heavy- ion collisions 0,^]. This has been incorporated into the BUU and QMD approaches 
based on various phenomenological parameterizations that are constrained by the empirical 
information from proton- nucleus scattering experiments [p|,P,p!0| . 

On the other hand, the relativistic transport models are usually based on the 

Walecka-type effective Lagrangians [|T1|]. In this model, nucleons feel both attractive scalar 
and repulsive vector potentials. The corresponding Schrodinger equivalent potential, that is 
to be compared with the potentials used in BUU and QMD calculations, is automatically 
and explicitly momentum dependent. In the mean field approximation both the scalar and 
vector potentials are momentum independent. The Schrodinger equivalent potential then 
increases linearly with nucleon kinetic energy, which is in disagreement with what has been 
learned from the Dirac phenomenology analysis of proton-nucleus scattering at intermediate 
energies 
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In a recent paper, we have shown that the explicit momentum dependence in nucleon 
scalar and vector potential is important for the description of nucleon flow in heavy-ion 
collisions at 1 GeV/nucleon region ||13|. In that work, the momentum dependence cor- 
rection was introduced by scaling the momentum independent mean field potential with a 
momentum dependent function extracted from the Dirac phenomenology. There have been 



similar attempts to reproduce the correct momentum dependence of the Schrodinger equiv- 
alent potential by including explicit momentum dependence in the nucleon scalar and vector 
potentials |T^. 

A main purpose of this paper is to study microscopically the origin of the momentum 



dependence in the scalar and vector potentials used in relativistic heavy-ion collisions [ITS 



The self-consitent Dirac Brueckner (DBHF) approach |l5Hl9| is employed. It is well-known 
that the DBHF approach reproduces nicely the nuclear matter properties starting from re- 
alistic nucleon-nucleon interactions that are constrained by the two-nucleon data. In most 
DBHF calculations the scalar and vector potentials are assumed to be momentum inde- 
pendent. This is a good approximation for equilibrium nucleon matter, where the nucleon 
momenta involved are not very large. The momentum dependence was studied by ter Haar 



and Malfliet |16], and they found a relatively weak momentum dependence for both the real 



and imaginary parts of the scalar and vector potentials. Recently, Sehn et al. studied 
the momentum dependence of the nucleon scalar and vector potentials in the first-order 
approximation. They first carried out the momentum independent DBHF calculation as 
in Ref. |T^. The scalar and vector potential so obtained were identified with a nucleon at 
Fermi surface, namely |k| = kp. The potentials for momenta below and above kp were 



then calculated using the formalism of Horowitz and Serot |^ . The self-consistency in the 
DBHF calculation has thus been neglected in Ref. [^]. We will compare our results to those 
of Refs. 0,0 . 



In this work, we follow essentially the DBHF formalism of Ref. [|T^. The basic quantity 
is the G matrix that satisfies the in-medium Thompson equation. 



d^k 

G{q',q\P,~z) = V{q',q)+ I J^,V{q\k) 



fm{k)V Q{k,P) 



\E{k)J 2E{q)-2E{k) 



G{k,q\P,2) 



where E{k) = \Jm? + (-P/2 + A;)^, P is cm. momentum, z = 2E{q) and m{k) = m + Us{k). 
The single-particle potential can be calculated once the G matrix is known. 



E,,(*) = Re I gg) gg) (k, \G(i)\ k, - ,k 



(2) 



This potential is usually parameterized in terms of scalar and vector potential, 

nk)^'^Us{k) + Uy{k). (3) 
E[k) 



In Ref. |]T^, both the scalar and vector potentials are assumed to be momentum independent. 



Thus for a given density, these potentials can be determined by choosing two momenta in 
Eq. (3), usually Q.lkp and kp. In the momentum dependent case as we consider here, 
we have two unknowns {Us and Uy) for each given momentum, which cannot be uniquely 
determined without introducing further approximation. The customary way to do this is 
to assume some functional dependence on momentum for these potentials. For example, in 
Ref. it was assumed that Usy ~ ^sy + Bsyk"^ ■ Since the real part of the single particle 
potential vanishes at very high momentum, we expect the following parameterization to be 
more appropriate 

1 + asyyk/kFr^-'^ 



Actually, in Ref. ||T3[, this kind of functional dependence was found to fit the Dirac phe- 
nomenology rather well. One of the purposes of this work is thus to study whether this kind 
of simple parameterization is supported by microscopic DBHF calculation. Numerically we 
find that it is rather difficult to determine the six parameters altogether in the minimization 
procedure. We have checked the variation of (3sy around 0.5, and find that (3sy ~ 0.5 pro- 
vides a smooth self-consistent procedure, and the final results are not very sensitive to some 
small change in (3sy- The remaining four parameters are then determined by a minimization 
procedure, 

Sk{j:DB{k)-nk)f . (5) 

Jo 

Our self-consistent procedure is the following. For a given density as specified by Fermi 
momentum kp, we choose reasonable initial values for {Usy)o and asy- These, together with 
Eq. (4) provide Usy for all momenta. These potentials are then used as inputs in Eq. (1) to 
determine the G matrix, from which we obtain the single-particle potential S^s using Eq. 



(2). Our new set of {Usy)o and asy can then be determined by minimizing Eq. (H) against 
TiDB- This procedure is continued until desired self- consistency is achieved. The quality of 
our self-consistency procedure can be seen in Fig. |^ where we compare the single-particle 
potential from the DBHF calculation (solid symbols) and from Eqs. (3) and (4). There are 
slight deviations near k = because in our definition of x^, Eq- (5), the contributions from 
small momenta are suppressed. The self consistent parameters are summarized in Tab. |. 

In Fig. 1^ we show our results for the momentum dependence of the scalar and vector 
potentials at normal nuclear matter density. Both are seen to decrease in magnitude as 
nucleon momentum increases. Also shown in the figure are the results of ter Haar and 
Malfiiet (HM) |[T^ and those of Sehn et al. (SFF) [|20|. Our results are very close to those 
of ter Haar and Malfiiet, which supports our parametrization functional, Eq. (^. The 
results of Sehn et al. show a much stronger momentum dependence than ours and those of 
ter Haar and Malfiiet |jl6[, because of the lack of the self-consistency. We show in Fig. |^ 



the momentum dependence of the nucleon scalar and vector potentials at several different 
densities. 

From the scalar and vector potential we can calculated the Schrodinger equivalent po- 
tential given by 

\ m J 2m 

where Ekin is the nucleon kinetic energy. The results are shown in Fig. |. The solid 
symbols are the Schrodinger equivalent potential extracted from the Dirac phenomenology 
analysis [|12[ . The solid and dashed curves are obtained with and without explicit momentum 
dependence in the scalar and vector potentials. The dotted is from ref . [|16[ . It is seen that 
the explicit momentum dependence in the scalar and vector potentials improves, to some 
extent, the agreement with the Dirac phenomenology. The Schrodinger equivalent potential 
from our self-consistent calculation still show stronger kinetic energy dependence than that 
of the Dirac phenomenology. The onset of inelastic processes at high beam energies might 
explain part of this discrepancy. 



It is also of interest to see how the nuclear equation of state are affected by the explicit 
momentum dependence. 

E/A ^ r f + '-m] - M. (7) 

' Jo^"- <Pk Jo V E{k) 2 ^ V 

In Fig.|, the binding energies are plotted as a function of density. The solid and dashed 
curves are obtained with and without explicit momentum dependence in the scalar and vector 
potentials. It is seen that the effects of the explicit momentum dependence on equilibriu 
nuclear matter is insignificant. The saturation density in both cases are about ~ 0.18 fm~^. 
The binding energy decreases by about 0.2 MeV when explicit momentum dependence is 
included. 

In summary, we studied the momentum dependence of nuclear scalar and vector poten- 
tials self-consistently in the formalism of Dirac Brueckner approach. The magnitudes of both 
potentials have been found to decrease slowly with increasing nucleon momentum or kinetic 
energy. This improves, to some extent, the agreement with the Dirac phenomenology. The 
effects of explicit momentum dependence on the properties of equilibrium nuclear matter 
were found to be small. 
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TABLES 



uip/po) 


UsoiMeV) 




U,o{MeV) 


ay 


0.70 


-426.2505 


0.6178 


370.6928 


0.8754 


0.85 


-474.3686 


0.4283 


408.2838 


0.5805 


1.00 


-494.2272 


0.3426 


420.5226 


0.4585 


1.12 


-499.1539 


0.2790 


420.0603 


0.3687 


1.24 


-518.8160 


0.2566 


436.9299 


0.3361 


1.37 


-537.0683 


0.2357 


453.5746 


0.3050 


1.51 


-586.0774 


0.2652 


508.8679 


0.3515 



TABLE I. Self-consistent parameters with po = 0.166 fm 



9 



FIGURES 



-20 - 



> 

CD 



CO 
•I— ( 

fi 

0) 

o 



-40 - 



-60 



-80 



u=1.51 / 



u-0.70 




^Xu=1.00 







0.5 



1 



1.5 



k (fm-i) 

FIG. 1. Momentum dependence of single particle potential. The solid symbols are the results 



of DBHF calculation and the lines are from Eqs. (3) and (4). 
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FIG. 2. nucleon scalar and vector potential as a function of nucleon momentum at normal 
nuclear matter density. The solid curves are our results, while the dashed and dotted curve are 
the results from Refs. pi (HM) and pO] (SFF), respectively. 
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FIG. 3. nucleon scalar and vector potential as a function of nucleon momentum at several 



different densities. 
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FIG. 4. Schrodinger equivalent potential in the DBHF calculation with and without explicit 
momentum dependence. The solid circles are extracted from the Dirac phenomenology analysis 
|ll. The results of Haar and Malfliet (HM) pi are plotted in dotted line. 
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FIG. 5. Binding energies in the DBHF calculation with and without explicit momentum de- 
pendence. 
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